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Abstract 

We extend a recent three-loop calculation of nuclear matter in chiral perturbation 

theory by including the effects from two-pion exchange with single and double virtual 

A(1232)-isobar excitation. Regularization dependent short-range contributions from pion- 

loops are encoded in a few NN-contact coupling constants. The empirical saturation point 

^^ ' of isospin-symmetric nuclear matter, £"0 = — 16MeV, po = 0.16 fm~'^, can be well re- 

f— ^ . produced by adjusting the strength of a two-body term linear in density (and tuning an 

CN i emerging three-body term quadratic in density). The nuclear matter compressibility comes 

out as K = 304 MeV. The real single-particle potential U{p, kjo) is substantially improved 

/^ ■ by the inclusion of the chiral vrA^A-dynamics: it grows now monotonically with the nu- 

cleon momentum p. The effective nucleon mass at the Fermi surface takes on a realistic 

value of M*{kfQ) = 0.88M. As a consequence of these features, the critical temperature 

p^^ ' of the liquid-gas phase transition gets lowered to the value Tc — 15 MeV. In this work we 

^ , continue the complex-valued single-particle potential U{p, kf) + i W{p, kf) into the region 

^f^ I above the Fermi surface p > kf. The effects of 27r-exchange with virtual A-excitation on 

^^ ■ the nuclear energy density functional are also investigated. The effective nucleon mass 

^^ . associated with the kinetic energy density is M*{pq) = 0.64M. Furthermore, we find that 

T-|- \ the isospin properties of nuclear matter get significantly improved by including the chi- 

CD ' ral vrA^A-dynamics. Instead of bending downward above pQ as in previous calculations, 

the energy per particle of pure neutron matter En{kn) and the asymmetry energy A{kf) 

now grow monotonically with density. In the density regime p = 2pn < 0.2 fm~^ rele- 

Q ' vant for conventional nuclear physics our results agree well with sophisticated many-body 

2 ■ calculations and (semi)-empirical values. 
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1 Introduction and preparation 

In recent years a novel approach to the nuclear matter problem based on effective field theory 
(in particular chiral perturbation theory) has emerged. Its key element is a separation of long- 
and short- distance dynamics and an ordering scheme in powers of small momenta. At nuclear 
matter saturation density po — 0.16 fm~^ the Fermi momentum kfo and the pion mass m^ 
are comparable scales (/c/o — 2mn), and therefore pions must be included as explicit degrees of 
freedom in the description of the nuclear many-body dynamics. The contributions to the energy 
per particle E{kf) of isospin-symmetric nuclear matter as they originate from chiral pion-nucleon 
dynamics have been computed up to three- loop order in refs.^|2]. Both calculations are able 
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to reproduce correctly the empirical saturation point of nuclear matter by adjusting one single 
parameter (either a coupling go+gi — 3.23 [T] or a cut-off A ~ 0.65 GeV j2]) related to unresolved 
short-distance dynamics. The novel mechanism for saturation in these approaches is a repulsive 
contribution to the energy per particle generated by Pauli-blocking in second order (iterated) 
one-pion exchange. As outlined in section 2.5 of ref.[2, this mechanism becomes particularly 
transparent by taking the chiral limit 771.,^ = 0. In that case the interaction contributions to 
the energy per particle are completely summarized by an attractive fcj-term and a repulsive 
ki-term where the parameter-free prediction for the coefficient of the latter is very close to the 
one extracted from a realistic nuclear matter equation of state. 

The single-particle properties, represented by a complex-valued momentum and density de- 
pendent nucleon selfenergy U{p,kf)+i W{p, kj), have been computed within our approach [2] in 
ref.jHJ. The resulting potential depth f/(0, kjo) = — 53.2MeV is in good agreement with that of 
the empirical nuclear shell 0] or optical model i5j. However, the momentum dependence of the 
real single-particle potential U{p, kjo) with its up- and downward bending (see Fig. 3 in ref.j2]) 
turns out to be too strong. As a consequence, the nominal value of the effective nucleon mass 
at the Fermi surface p = kfo would be much too high: M*{kfo) ~ 3M. On the other hand, the 
single-particle properties around the Fermi surface are decisive for the spectrum of thermal exci- 
tations and therefore they crucially influence the low temperature behavior of isospin-symmetric 
nuclear matter. The rather high critical temperature T^ ~ 25.5 MeV for the liquid-gas phase 
transition obtained in ref. [Hj is a visible manifestation of this intimate relationship. 

While there is obviously room and need for improvement in our approach, one must also 
note at the same time that the single-particle properties in the scheme of Lutz et al. IJ (where 
explicit short-range terms are iterated with pion-exchange) come out completely unrealistic [7j . 
The potential depth of ^(0, fc/o) = ~20 MeV is by far too weakly attractive. Most seriously, 
the total single-particle energy Tkin(p) + U{p, /c/o) does not rise monotonically with the nucleon 
momentum p, thus implying a negative effective nucleon mass at the Fermi surface (see Fig. 3 
in ref.|[73). This ruins the behavior of nuclear matter at finite temperatures since a critical 
temperature of Tc > 40 MeV exceeds acceptable values by more than a factor of two. 

The isospin properties of nuclear matter have also been investigated in ref. [2j. The prediction 
for the asymmetry energy at saturation density A{kfo) = 33.8 MeV is in good agreement with 
its empirical value. However, one finds a downward bending of A{kf) at densities p > 0.2 fm~^. 
Such a behavior of the asymmetry energy A{kf) is presumably not realistic. The energy per 
particle of pure neutron matter En{kn) as a function of the neutron density p„ = kf^/Sii^ shows 
a similar downward bending behavior (see Fig. 8 in ref.pi) and at lower neutron densities, there 
is only rough agreement with realistic neutron matter calculations. The mere fact that neutron 
matter came out to be unbound in ref. [2j with no further adjusted parameter was however non- 
trivial. The isospin properties of nuclear matter in the scheme of Lutz et al. jT] with a second 
free parameter adjusted, are qualitatively the same. The dashed curves in Figs. 6,7 of ref. [Zj 
display a downward bending of En{kn) and A{kf) at even lower densities, p > 0.15 fm~^. An 
extended version of that approach with pion-exchange plus two zero-range NN-contact interac- 
tions iterated to second order and in total four adjustable parameters has been studied recently 
in ref. [H| . The finding of that work is that within such a complete fourth order calculation (thus 
exhausting all possible terms up-to-and-including 0{kj)) there is no optimal set of the four 
short-range parameters with which one could reproduce simultaneously and accurately all semi- 
empirical properties of nuclear matter. The conditions for a good neutron matter equation of 
state and equally good single-particle properties (and consequently a realistic finite temperature 
behavior) are in fact mutually exclusive in that approach. 

Calculations of nuclear matter based on the universal low-momentum nucleon-nucleon po- 
tential ^ow-k have recently been performed in ref. jH] . The results obtained so far in Hartree-Fock 



or Brueckner-Hartree-Fock approximation are unsatisfactory (see Figs. 1,2 in ref. [Oj), since no 
saturation occurs in the equation of state. It has been concluded that for the potential Viow-k the 
Brueckner-Hartree-Fock approximation is applicable only at very low densities. These findings 
together with the identification of a successful saturation mechanism in the chiral approaches 
hint at the fact that the Brueckner ladder does not generate all relevant medium modifications 
which set in already at very low densities (typically at about one-tenth the equilibrium density 
of nuclear matter, corresponding to Fermi momenta around kj ^ wi^). 

Up to this point the situation can be summarized as follows. Chiral two-pion exchange 
restricted to nucleon intermediate states (basically the second-order spin-spin and tensor force 
plus Pauli blocking effects), together with a single (finetuned) contact-term representing short- 
distance dynamics, is already surprisingly successful in binding and saturating nuclear matter 
and reproducing key properties such as the asymmetry energy and the compression modulus. 
However, the detailed behavior of the nucleon single-particle potential and the density of states 
at the Fermi surface are not well described at this order in the small-momentum expansion. 

Let us also comment on the relationship between our appoach to nuclear matter and the 
effective field theory treatment of (free) NN-scattering [HIl HH IE] • We are including to the 
respective order the same long-range components from (one- and) two-pion exchange. The short- 
range NN-contact terms are however treated differently. Instead of introducing these terms in a 
potential which is then iterated in a Lippmann-Schwinger equation, we adjust their strengths (in 
a perturbative calculation) to a few empirical nuclear matter properties. The resulting values 
of the coupling constants are therefore not directly comparable with those of the NN-potential. 
We note also that the important role of chiral two-pion exchange for peripheral NN-scattering 
is well established HH [HI CS] • 

From the point of view of the driving pion-nucleon dynamics the previously mentioned chiral 
calculations of nuclear matter [HllllHj are indeed still incomplete. They include only (S- and) P- 
wave Born terms but leave out the excitation of the spin-isospin-3/2 A(1232)-resonance, which 
is the prominent feature of low-energy vrA^-scattering. It is also well known that the two-pion 
exchange between nucleons with excitation of virtual A-isobars generates most of the needed 
isoscalar central NN-attraction. In phenomenological one-boson exchange models this part of 
the NN-interaction is often simulated by a fictitious "a" -meson exchange. A parameter-free 
calculation of the isoscalar central potential Vc{r) generated by 27r-exchange with single and 
double A-excitation in ref . [,.2j (see Fig. 2 therein) agrees almost perfectly with the phenomeno- 
logical " cr" -exchange potential at distances r > 2fm, but not at shorter distances. The more 
detailed behavior of the 27r-exchange isoscalar central potential with single virtual A-excitation 
is reminiscent of the van-der-Waals potential. It has the form 



Q 4 —2x 

with X = rriT^r and the prefactor includes the spin-isospin (axial) polarizability of the nucleon 
[T3], gl/f^A = 5.2 fm^ from the virtual A^ -^ A(1232) -^ N transition. The familiar r"^- 
dependence of the non-relativistic van-der-Waals interaction emerges in the chiral limit, m.,^ = 0. 
A consideration of mass scales also suggests to include the A(1232)-isobar as an explicit de- 
gree of freedom in nuclear matter calculations. The delta-nucleon mass splitting of A = 293 MeV 
is comparable to the Fermi momentum kfo ^ 262 MeV at nuclear matter saturation density. 
Propagation effects of virtual A(1232)-isobars can therefore be resolved at the densities of in- 
terest. Based on these scale arguments we adopt a calculational scheme in which we count the 
Fermi momentum kf, the pion mass 771^^ and the AA^-mass splitting A simultaneously as "small 
scales" . The non-relativistic treatment of the nuclear matter many-body problem naturally goes 



conform with such an expansion in powers of small momenta. Relativistic corrections are rele- 
gated to higher orders in this expansion scheme. The leading contributions from 27r-exchange 
with virtual A-excitation to the energy per particle (or the single-particle potential) are gener- 
ically of fifth power in the small momenta {kf,mT^,A). With respect to the counting in small 
momenta the effects from irreducible 27r-exchange evaluated in ref. |2l El E] belong to the same 
order. However, since the vrA^A-coupling constant is about twice as large as the vrA^A^-coupling 
constant one can expect that the A-driven 27r-exchange effects are the dominant ones. The im- 
portance of A(1232)-degrees of freedom has also been pointed out in the "ab-initio" calculations 
of the Illinois group [121 HI] . 

The purpose of the present paper is to present a calculation of (isospin-symmetric and 
isospin-asymmetric) nuclear matter which includes systematically all effects from 27r-exchange 
with virtual A-excitation up to three-loop order in the energy density. The contributions to 
the energy per particle (or the single-particle potential) can be classified as two-body terms 
and three-body terms. Two-body terms can be directly expressed through the NN-scattering 
T- matrix (i.e. the NN-potential in momentum space). Three-body terms on the other hand can 
be interpreted as Pauli-blocking effects on the two-body terms imposed by the filled Fermi-sea 
of nucleons. The notion of "three-body term" is taken here in a more general context, namely 
in the sense that three nucleons in the Fermi sea participate in interactions. The NN T-matrix 
generated by the in general ultra-violet divergent pion-loop diagrams requires regularization 
(and renormalization) . We adopt here a suitably subtracted dispersion-relation representation 
of the T-matrix where this procedure is accounted for by a few subtraction constants. The latter 
constants are understood to encode unresolved short-distance NN-dynamics. The associated k^- 
and /cj-terms in the energy per particle are then adjusted to some empirical property of nuclear 
matter (e.g. the maximal binding energy of 16MeV). 

Our paper is organized as follows: In section 2 we start with the equation of state of isospin- 
symmetric nuclear matter. We evaluate analytically the three-loop in-medium diagrams gener- 
ated by the chiral vriVA-dynaniics and perform the necessary adjustment of short-range parame- 
ters. Section 3 deals with the real single-particle potential U{p, kf) whose improved momentum 
dependence turns out to be a true prediction. In section 4 we reconsider the imaginary single- 
particle potential W{p, kf) of ref.0, but now extended into the region above the Fermi surface 
p > kf. Section 5 is devoted to the effective nucleon mass M*{p) and the strength function of 
the (Vp)^-term in the nuclear energy density functional. In section 6 we extend our calculation 
of isospin-symmetric nuclear matter to finite temperatures T. The main interest lies there on 
the critical temperature T^ of the first-order liquid-gas phase transition for which we find an 
improved value of Tc ~ 15 MeV. Sections 7, 8 and 9 deal with the equation of state of pure 
neutron matter En{kn), the asymmetry energy A{kf) and the isovector single-particle potential 
Ui{p, kf). These three quantities reveal the isospin properties of the underlying TriVA-dynamics. 
Explicit inclusion of the A(1232)-degree of freedom leads to a substantial improvement: the no- 
torious downward bending of E^ikn) and A{kf) observed in previous chiral calculations is now 
eliminated. Finally, section 10 ends with a summary and some concluding remarks. 

2 Equation of state of isospin-symmetric nuclear matter 

We start the discussion with the equation of state of isospin-symmetric nuclear matter for 
which one has a fairly good knowledge of the empirical saturation point. We first write down 
the contributions to the energy per particle E{kf) as they arise from 27r-exchange with single 
and double virtual A-isobar excitation. 

Fig. 1 shows the relevant one-loop triangle, box and crossed box diagrams contributing to 
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Figure 1: One-loop two-pion exchange diagrams with single and double A(1232)-isobar excita- 
tion. Diagrams for which the role of both nucleons is interchanged are not shown. 



the NN T-matrix (or the momentum space potential). The finite parts of these diagrams have 
been evaluated analytically in section 3 of ref . |E] employing the usual non-relativistic A *^ ttN 
transition vertices and A-propagator (see eq.(4) in ref.jHj). By closing the two open nucleon 
lines to either two rings or one ring one gets (in diagrammatic representation) the Hartree or Fock 
contribution to the energy density of nuclear matter. The Hartree contribution to the energy 
per particle evidently goes linear with the nucleon density p = 2fc^/37r^, namely E{kf)^'^^^ = 



— Vc'(O) p/2 with Vc'(O) the isoscalar central NN-amplitude at zero momentum transfer ^^. The 
Fock contribution on the other hand is obtained by integrating the spin- and isospin-contracted 
T-matrix (depending on the momentum transfer variable \pi — P2I) over the product of two 
Fermi spheres \'Pi^2\ < kf of radius kf. We separate regularization dependent short-range parts 
in the T-matrix (originating from the divergences of the loop diagrams) from the unique long- 
range terms with the help of a twice-subtracted dispersion relation. The resulting subtraction 
constants give rise to a contribution to the energy per particle of the form: 
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where -B3 and B^ are chosen for convenience as dimensionless. M = 939 MeV stands for the 
(average) nucleon mass. Note that eq.(l) is completely equivalent to the contribution of a mo- 
mentum independent and p^-dependent NN-contact interaction. We interpret the parameters 
i?3,5 to subsume all unresolved short-distance NN-dynamics relevant for isospin-symmetric nu- 
clear matter at low and moderate densities. The long-range parts of the 27r-exchange (two-body) 
Fock diagrams can be expressed as: 
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where ImVc, ImVFc, ImVr and ImPVV are the spectral functions of the isoscalar and isovector 
central and tensor NN-amplitudes, respectively. Explicit expressions of these imaginary parts 
for the contributions of the triangle diagram with single A-excitation and the box diagrams with 
single and double A-excitation can be easily constructed from the analytical formulas given in 
section 3 of ref. [131. The p- and /c/-dependent weighting function in eq.(2) takes care that 
at low and moderate densities this spectral integral is dominated by low invariant Trvr-niasses 
2r7i7r < p < 1 GeV. The contributions to the energy per particle from irreducible 27r-exchange 
(with only nucleon intermediate states) can also be cast into the form eq.(2). The corresponding 







Figure 2: Hartree and Fock three-body diagrams related to 27r-exchange with single virtual 
A-isobar excitation. They represent interactions between three nucleons in the Fermi sea. In 
the case of isospin-symmetric nuclear matter the isospin factors of these diagrams are 8, 0, and 
8, in the order shown. The combinatoric factor is 1 for each diagram. 



non- vanishing spectral functions read J13j : 
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The dispersion integrals J^^ dfilYa{. . .) in this and all following sections are understood to 
include the contributions from irreducible 27r-exchange (with only nucleon intermediate states). 
Next, we come to the three-body terms which arise from Pauli blocking of intermediate 
nucleon states (i.e., from the —2TTd{kf — \p\) terms of the in-medium nucleon propagators ^). 
The corresponding closed Hartree and Fock diagrams with single virtual A-excitation are shown 
in Fig. 2. In the case of isospin-symmetric nuclear matter their isospin factors are 8,0, and 8, 
in the order shown. For the three-loop Hartree diagram the occurring integral over the product 
of three Fermi spheres of radius kf can be solved in closed form and the contribution to the 
energy per particle reads: 
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with the abbreviation u = kf/rriTt where m^r = 135 MeV stands for the (neutral) pion mass. 
The delta-propagator shows up in this expression merely via the (reciprocal) mass-splitting 
A = 293 MeV. Additional corrections to the delta-propagator coming from differences of nu- 
cleon kinetic energies etc. will make a contribution at least one order higher in the small- 
momentum expansion. In eq.(5) we have already inserted the empirically well-satisfied relation 
5'7rAfA = 3(77rAr/"\/2 for the vrA^A-coupling constant together with the Goldberger-Treiman re- 
lation ^f^Tv = Qa^I f-K (see e.g. eq.(5) in ref. 114. for the A — >• Nn decay width). As usual 
/tt = 92.4 MeV denotes the weak pion decay constant and we choose the value qa = 1-3 in order 
to have a pion-nucleon coupling constant of g.„M = 13.2 [TH]. Via the parameter C, we have 
already included in eq.(5) the contribution of an additional three-nucleon contact interaction 
~ {(g\/Af^) [NN)^. One notices that the 27r-exchange three-body Hartree diagram shows in 
the chiral limit m^r = the same quadratic dependence on the nucleon density. In that limit the 
momentum dependent vrA^A-interaction vertices get canceled by the pion-propagators and thus 
one is effectively dealing with a zero-range three-nucleon contact-interaction. It is important to 
point out that this equivalence holds only after taking the spin-traces but not at the level of the 
(spin- and momentum dependent) 27r-exchange three-nucleon interaction. The contribution of 
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the (right) three-body Fock diagram in Fig. 2 to the energy per particle reads: 
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Evidently, the three-body Fock term in eq.(6) is attractive. The remaining contributions to 
E[kf) from the (relativistically improved) kinetic energy, from the Ivr-exchange Fock diagram 
and from the iterated l7r-exchange Hartree and Fock diagrams have been written down in 
eqs.(5-ll) of ref.j2j. The strongly attractive contribution from iterated Ivr-exchange linear in 
the density and the cutoff A (see eq.(15) in ref.pj) is now of course not counted extra since 
B^ in eq.(l) collects all such possible terms. Adding all pieces we arrive at the full energy per 
particle E{kf) at three-loop order. It involves three parameters, B^ and B^ of the two-body 
contact term eq.(l) and ( which controls a three-body contact term in eq.(5). 

Let us first look at generic properties of the nuclear matter equation of state in our calcu- 
lation. Binding and saturation occurs in a wide range of the two adjustable parameters -83,5. 
However, with the full strength {( = 0) of the repulsive p^-term from the 27r-exchange three- 
body Hartree diagram (see eq.(5)) the saturation curve rises much too steeply with increasing 
density. This causes a too low saturation density po and a too high nuclear matter compress- 
ibility, K > 350 MeV. We cure this problem in a minimal way by introducing an attractive 
three-body contact term. With ( = —3/4 the remaining repulsive p^-term in eq.(5) gets can- 
celed by an analogous attractive contribution from the three-body Fock diagram. Clearly, the 
need for introducing an attractive three-body contact term into our calculation points to some 
short- distance physics whoose dynamical origin lies outside the present framework of pertur- 
bative chiral pion-nucleon interactions. It will become clear in the following sections that the 
predictive power of our calculation is nevertheless not reduced by this procedure. 

We fix the minimum of the saturation curve E{kf) to the value ^0 = —16.0 MeV. With B^ 
adjusted to the value -B3 = —7.99 and B^ taken to be zero, B^ = 0, the full line in Fig. 3 results. 
The predicted value of the saturation density is po = 0.157 fm~^, corresponding to a Fermi 
momentum of fc/o = 261.6 MeV = 1.326 fm~^. This is very close to the semi-empirical value 
Po = 0. 158 fm~^ obtained by extrapolation from inelastic electron scattering off heavy nuclei ^H] ■ 
The decomposition of the negative binding energy Eq = —16.0 MeV into contributions of second, 
third, fourth and fifth power in small momenta reads: Eq = (21.9 — 145.5 + 107.8 — 0.2) MeV with 
the typical balance between large third and fourth order terms j2]. The very small fifth order 
term splits furthermore as (—13.8 + 13.6) MeV into the contribution from the three-body contact- 
interaction (proportional to C = —3/4) and a remainder. Evidently, since Eq = —16.0 MeV 
is a small number that needs to be finetuned in our calculation there remains the question 
of the "convergence" of the small momentum expansion. The nuclear matter compressibility 
K = kj-^E'^kfo) related to the curvature of the saturation curve at its minimum comes out 
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Figure 3: The energy per particle E{kf) of isospin-symmetric nuclear matter as a function of 
the nucleon density p = 2/cj/37r^. The dashed line refers to the result of ref.p], with only pions 
and nucleons as active degrees of freedom. The full line includes effects from 27r-exchange with 
virtual A-excitation. The short-range parameters are B3 = —7.99 and -B5 = 0. 



as K = 304 MeV. This number is somewhat high but still acceptable, since it exceeds e.g. 
the value K = 272 MeV obtained in the relativistic mean-field model of ref. [201 only by 12%. 
Inspection of Fig. 3 shows that the saturation curve is well approximated by a shifted parabola 
E{kf) = Eq{2pq — p)p/pI with a second zero-crossing at twice nuclear matter density 2pQ. This 
leads to a compressibility estimate of i^ = — 18^o = 288 MeV, not far from the calculated value. 
We have also studied variations of the parameter B^. Within the limits set by a stable saturation 
point (and realistic single-particle properties, see next section), the effects on the compressibility 
K are marginal (less than 10 MeV reduction). Therefore we stay with the minimal choice -B5 = 
(together with B-^ = —7.99 and ( = —3/4). The dashed line in Fig. 3 shows for comparison 
the equation of state resulting from our previous chiral calculation f2^ with no vrA^A-dynamics 
included. In that work the saturation density po = 0.178 fm~^ came out somewhat too high, 
but the compressibility K = 255 MeV had a better value. The stronger rise of the full curve in 
Fig. 3 with density p is a consequence of including higher order terms in the (small-momentum) 
fcj-expansion. 



3 Real part of single-particle potential 

In this section we discuss the real part U{p, kf) of the single-particle potential. As outlined in 
ref.|H] the contributions to the (real) nuclear mean-field U{p,kf) can be classified as two-body 
and three-body potentials. The parameters -83^5 introduced in eq.(l) reappear in a contribution 
to the two-body potential of the form: 



U2{p,kj 



\ict) 



2B. 



l3 

M2 



+ B, 



7,3 

^/ 

3M4 



{3kj + 5p^) 



(9) 



Its density- and momentum-dependence is completely fixed by the Hugenholtz- Van-Hove theo- 
rem j2Ij and a sum rule which connects it to the energy per particle E{kf)^'^^ (see eqs.(5,7) in 
ref.jH]). The Fock diagrams of 27r-exchange with virtual A-excitation give rise to a contribution 
to the two-body potential which can be written as a (subtracted) dispersion integral: 
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By opening a nucleon line in the three-body diagrams of Fig. 2 one gets (per diagram) three 
different contributions to the three-body potential. In the case of the (left) Hartree diagram 
they read altogether: 
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with the abbreviation x = p/m.^. Note that the paramter ( = —3/4 related to the three- 
body contact-interaction has no influence on the momentum dependence of the single-particle 
potential. On the other hand the (right) Fock diagram in Fig. 2 generates a total contribution 
to the three-body potential of the form: 



U3{p,kjY^^ 



9>U ' i2Gl(x,u) + Gl(x,u) 



4A(47r^ 



•^SK-^^ 



+ / d^ 

'0 



465(4, u) 7; h 2Gt(4, u) 



dx 



dx 



fl2) 



with Gs,t{x,u) defined in eqs.(7,8). The real single-particle potential U{p,kf) is completed by 
adding to the terms eqs.(9-12) the contributions from Ivr-exchange and iterated Ivr-exchange 
written down in eqs.(8-13) of ref.jB]. The slope of the real single-particle potential U{p, kf) at 
the Fermi surface p = kf determines the effective nucleon mass (in the nomenclature of ref . j22] , 
the product of " fc-mass" and " E'-mass" divided by the free nucleon mass M = 939 MeV) via a 
relation: 

M dU{p, kf) 



M*{kf) = M 



2M2 ^ k 



f 



dp 



-i -1 



p=kf 



(13) 



The second term —k'j/2M'^ in the square brackets stems from the relativistic correction —p^/8M^ 
to the kinetic energy. 

The full line in Fig. 4 shows the real part of the single-particle potential U{p, kfo) at sat- 
uration density kfQ = 261.6 MeV as a function of the nucleon momentum p. The dotted line 
includes in addition the relativistically improved kinetic energy Tkin(p) = p^ /2M — p^/8M^. 
With the parameters i?5 = 0, C = —3/4 fixed and B^ = —7.99 adjusted in section 1 to the 
binding energy at equilibrium we find a potential depth of f/(0, kfo) = —78.2 MeV. This is very 
close to the result f/(0, fc/o) — —80 MeV of the relativistic Dirac-Brueckner approach of ref. 
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Figure 4: Full line: real part of the single-particle potential U{p, kfo) at saturation density 
kfQ = 261.6 MeV as a function of the nucleon momentum p. The dotted line includes in addition 



the relativistically improved kinetic energy T\^in{p) 
into the region above the Fermi surface p > kjQ. 



p /2M—p /8M . Both curves are extended 



For comparison, the calculation of ref. [23] based on the phenomenological Paris NN-potential 
finds a somewhat shallower potential depth of U{0, kfo) ^ —64 MeV. One observes that with the 
chiral TrA^A-dynamics included, the real single-particle potential U{p, fc/o) grows monotonically 
with the nucleon momentum p. The downward bending above p = 180 MeV displayed in Fig. 3 
of ref. |2| is now eliminated. The slope at the Fermi surface p = kfo translates into an effective 
nucleon mass of M*{kfQ) = 0.88M. This is now a realistic value compared to M*{kfQ) ~ 3M 
obtained in our previous calculation [HJ E] without any explicit A-isobars. Note also that the 
chiral approach of ref. [Hj (where both explicit short-range terms and pion-exchange are iterated) 
has found the lower bound M*{kfQ) > 1AM. 

The dotted curve in Fig. 4 for the total single-particle energy Tkin(p) + f/(p, fc/o) hits the value 
E{kfQ) = Eq = —16 MeV at the Fermi surface p = kfQ, as required by the Hugenholtz- Van-Hove 
theorem |2l]. This important theorem holds strictly in our (perturbative) calculation, whereas 
(non-perturbative) Brueckner-Hartree-Fock approaches often fail to respect it |2S1- In Fig. 4 we 
have also extended both curves into the region above the Fermi surface p > fc/o- In general 
this extension is not just an analytical continuation of the potential from below the Fermi 
surface. Whereas eqs.(9-12) for the contributions from 27r-exchange with A-excitation apply in 
both regions, there are non-trivial changes in the expressions from iterated l7r-exchange. These 
modifications are summarized in the Appendix. The smooth rise of U{p, kfo) as it crosses the 
Fermi surface and proceeds up to p ~ 400 MeV is compatible with other calculations [23 121] • 
Beyond this momentum scale one presumably exceeds the limits of validity of the present chiral 
perturbation theory calculation of nuclear matter. 

The full line in Fig. 5 shows the potential depth f/(0, kf) for a nucleon at the bottom of 
the Fermi sea as a function of the nucleon density p = 2/cj/37r^. The band spanned by the 
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Figure 5: The full line shows the real part of the single-particle potential f/(0, kf) at nucleon 
momentum p = versus the density p = 2kj/37T'^. The band is obtained from the universal 
low-momentum NN-potential l1ow-k in linear density approximation. 



dotted lines stems from the universal low-momentum NN-potential l^ow-k 
approximation. In this approximation the potential depth simply reads: 



f/(0, k 



f) 



^CSo) 



Svrp 
2M 



K, 



i'So) 



low—k 



(o,o) + i^S(o,o) 



[2^ in linear density 



(14) 



with V^^'_'^^{0,0) ~ -1.9 fm and l/i^^iUO, 0) ~ -(2.2 ± 0.3) fm [201 EH EH], the two S-wave 
potentials at zero momentum. It is interesting to observe that both potential depths agree 
fairly well at low densities, p < 0.07 fm~^. This agreement is by no means trivial since l^ow-k is 
constructed to reproduce accurately the low-energy NN-scattering data (phase-shifts and mixing 
angles) whereas our adjustment of the i?3-term (linear in density) is made at saturation density 
Po = 0.157 fm~^. It is evident from Fig. 5 that a linear extrapolation does not work from zero 
density up to nuclear matter saturation density. Strong curvature effects set in already at Fermi 
momenta around kf ~ r7i,r if (and only if) pion-dynamics is treated explicitly.^ We note also 
that an "improved" determination of the potential depth f/(0, kj) from V^ow-k where one takes 
into account its momentum dependence in the repulsive Fock contribution leads to concave 
curves which bend below the straight dotted lines in Fig. 5. In the case V^iow-k(O'O) ~ — 1.9fm 
the potential depth ^(0, fc/o) at saturation density would increase to —132.6 MeV (compared to 
—113.4 MeV in linear density approximation). The present observations concerning the potential 
depth f/(0, kf) may indicate why calculations based on V^ow-k did so far not find saturation of 
nuclear matter jU]. It seems that the Brueckner ladder does not generate all relevant medium 
modifications which set in already at rather low densities kf ~ mT^ (if the pion-dynamics is 
treated explicitly). 

^As example for the extreme inherent non-linearities, consider the formula for the three-body potential 
1/3(0, kf)^^^ in eq.(ll). Its mathematical Taylor-series expansion converges only for kf < m,r/2. This corre- 
sponds to tiny densities, p < 0.0027 fm~^. 
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4 Imaginary part of single-particle potential 



In this section, we reconsider the imaginary part W{p, kf) of the single-particle potential. To 
the three-loop order we are working here it is still given completely by iterated Ivr-exchange 
with no contribution from the vrA^A-dynamics. The new aspect here is the extension into the 
region above the Fermi surface p > kf, which is not an analytical continuation from below the 
Fermi surface. As outhned in ref.j2| the contributions to W{p, kf) can be classified as two-body, 
three-body and four-body terms. From the Hartree diagram of iterated l7r-exchange one finds 
altogether (in the region above the Fermi surface p > kf): 



Wip,kf)^^'^ 



■ng\Mm l 
(4vrA 

+ 

V X 



f 9 + 6m^ + 2xM ln[l + {u + xf 



+ - 



X 



2x^-9- Qu^\ ln[l + (x - uf] + 4mx(2 - u^ 
(7 + 15m^ — 15x^)[arctan('U + x) — arctan(x — u)] 



I2u^ 21m 



-8M^ln(l + 4u2 
'2s^ + s^ 



15u' 



arctan 2u 



5 2 

fu/x 



xA, 



1 + S2 



-21n(l + s' 



+ 



dy Ai 



i"a/x 

+ 3e{V2u -x) dy{x^y^ - u^) 

3s^x\y^ - 1) 



1 + S2 



9s^ 



+10x1/(3 arctan s - 3s + s^) + (9 + 6m^ - Qx^y^)[s^ - ln(l + s^)] 



(15) 



with X = p/m.^ and the auxiliary functions y^ 



1 — M^/x^ and s = xy + ^/' 



u^ 



x^ + x^y^. 

In order to keep the notation compact we have introduced the antisymmetrization prescription 
•^y[f{y)] = fiy) ~ fi~y)- Note that there is a term in eq.(15) which vanishes identically above 
p = \/2kf. A geometrical explanation for this non-smooth behavior is that an orthogonal pair 
of vectors connecting the origin with two points inside a sphere ceases to exist if the center 
of the sphere is displaced too far from the origin (namely by more than \/2 times the sphere 
radius). The orthogonality of the (momentum difference) vectors is imposed here by the non- 
relativistic on-mass-shell condition for a nucleon. The combined two-body, three-body and 
four-body contributions to W{p, kf) from the iterated Ivr-exchange Fock diagram read on the 
other hand (for p > kf): 



W{p,kf)^^^ 



Txg\Mml 
(47r/.)4 



u 



{x+u)/2 



M"x + — + — / di (2^ - x) 

5X 2x J(^x-u)/2 ^ 

"1 Z)/i „,2 



U 



i+Ae 

1 + 2^2 



ln(l + 4^^ 



xA 



ymin 

e - in(i + e) 



+A0(72m-x) f dy r dz ^^ ^ ^^ ' =^A, 

■'" ■ -'J/min 



\ — y^ — Z^ 



s^ - ln(l + s^ 



ljjyj\e[Hi.^^)-^%i-{)], 



(16) 



with some new auxiliary functions t 



xz 



+ y/v? — x^ + x'^z'^ and <J = ^y + \/u^^-^^r^hp' 



and R = ^y{l + x^ — ^^)^ + 4:^^(1 — y"^). It is also interesting to consider the imaginary single- 
particle potential W{p, kf) in the chiral limit itIt^ = 0. One finds the following closed form 
expressions: 

,5 _ 

p>kf, (17) 



wip, kf)u., = g^{ 5 - kjp - A(2^J - /)V^ eiV2kf - p)] 
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Figure 6: The imaginary part of the single-particle potential W{p, kfo) at saturation density 
kfQ = 261.6 MeV as a function of the nucleon momentum p. The quadratic behavior around the 
Fermi surface p = kfo with a sign change of the curvature is required by Luttinger's theorem 



W{p,kf)l 



97rg\M . 2n2 ^ , 

[kr-p) , p<k 



/' 



4(47r^)4 

to which the iterated Ivr-exchange Hartree and Fock diagrams have contributed in the ratio 
4 : —1. The analytical results in eqs. (17,18) agree with Galitskii's calculation [2^ of a contact- 
interaction to second order. In the chiral limit m^ = the spin-averaged product of two ttNN- 
interaction vertices gets canceled by the pion propagators and thus one is effectively dealing with 
a zero-range NN-contact interaction at second order. The agreement with Galitskii's result [29] 
serves as an important check on the technically involved calculation behind eqs. (15,16). 

Fig. 6 shows the imaginary part of the single-particle potential W{p, kfo) at saturation den- 
sity kfo = 261.6 MeV as a function of the nucleon momentum p. The quantity ±2W{p,kf) 
determines the width of a hole-state or a particle-state of momentum p < kf or p > kf. The 
finite life time of these states originates from redistributing energy into additional particle-hole 
excitations. Our predicted value W^ (0, fc/o) = 24.0 MeV at p = lies in between the results 
W{0, kfo) ~ 20 MeV of ref.jHI] employing the Gogny Dl effective interaction and W{0, kfo) ~ 
40 MeV of ref. [21] using the Paris NN-potential. As a consequence of the decreasing phase- 
space for redistributing the hole-state energy, the curve in Fig. 6 drops with momentum p and 
W{p, kfo) reaches zero at the Fermi surface p = kfQ. According to Luttinger's theorem [HO] this 
vanishing is of quadratic order r^ {p — kf)"^, a feature which is clearly exhibited by the curve in 
Fig. 6. 

When crossing the Fermi surface the curvature of the imaginary single-particle potential 
W{p, kf) flips the sign. From there on a rapid fall to negative values sets in. In fact the width 
Tgp = —2W{p, kf) represents the spreading of a single-particle state above the Fermi surface into 
two-particle-one-hole states with growing phase space as p — kf increases. The range of validity 
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of the present chiral perturbation theory calculation is again expected to be p < 400 MeV. 
The rapid growth of T^p beyond such limited momentum scales is not shared by the results of 
refs.|2aEIl- 



5 Nuclear energy density functional 

The energy density functional is a general starting point for (non-relativistic) nuclear structure 
calculations within the framework of the self-consistent mean-field approximation [^^ . In this 
context effective Skyrme forces [SHI IHH IHSj have gained much popularity because of their analyt- 
ical simplicity and their ability to reproduce nuclear properties over the whole periodic table. In 
a recent work jSE] we have calculated the nuclear energy density functional which emerges from 
(leading and next-to-leading order) chiral pion-nucleon dynamics. The calculation in ref. [5E] 
included (only) the Ivr-exchange Fock-diagram and the iterated Ivr-exchange Hartree and Fock 
diagrams. These few components alone already lead to a good nuclear matter equation of 
state E{kf). Therefore the interest here is on the additional effects from 27r-exchange with 
virtual A-excitation contributing one order higher in the small momentum expansion. Going 
up to quadratic order in spatial gradients (i.e. deviations from homogeneity) the energy density 
functional relevant for N=Z even-even nuclei reads 



S[p,T]=pE{kf) + 



T 



P^ 



1 5fc/ 

2M 56M3 



+ F,{kf) 



+ (Vp)2Fv(M 



(19) 



with p{f) = 2k^r{f)/37c'^ the local nucleon density (expressed in terms of a local Fermi momen- 
tum kf{f)) and T{f) the local kinetic energy density. We have left out in eq.(19) the spin-orbit 
coupling term since the corresponding results (for A-driven 27r-exchange three-body spin-orbit 
forces^) can be found in ref. |SZ]. In phenomenological Skyrme parameterizations, the strength 
function Fr{kf) goes linearly with the density p, while F^{kf) is constant. The starting point 
for the construction of an explicit nuclear energy density functional £[p, r] is the bilocal density- 
matrix as given by a sum over the occupied energy eigenf unctions: Ylia&cc ^ai^ — a/2)\E'Jj(r + 
a/2). According to Negele and Vautherin [2H] it can be expanded in relative and center-of-mass 
coordinate, a and r, with expansion coefficients determined by purely local quantities (nucleon 
density, kinetic energy density and spin-orbit density). As outlined in section 2 of ref. ^j the 
Fourier transform of the (so-expanded) density matrix defines in momentum space a "medium 
insertion" for the inhomogeneous many-nucleon system: 



r(p,g; 



d^re-''^-''e{kf-\p\){ 1 + 



l^(5p--3*5) 



\pk] - ivv 



(20) 



The strength function F^lkf) in eq.(19) emerges via a perturbation on top of the density of 
states 9{kf — |p |). As a consequence of that, F^ikf) can be directly expressed in terms of the 
real single-particle potential U{p, kf) as: 



^^(*'' - 1 



dpp^ihp^ 



3kj) 



U{p,kf 



(21) 



Note that any p-independent contribution, in particular the i?3-term in eq.(9) and the (^-term 



in eq.(ll), drops out. In the medium insertion eq.(20) 



r 



3pkj/5 is accompanied by — V^p/4. 



''interestingly, this three-body spin-orbit couphng is not a relativistic effect but independent of the nucleon 
mass M . 
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After performing a partial integration one is lead to the following decomposition: 

where Ffi{kf) comprises all those contributions for which the (Vp)^-factor originates directly 
from the momentum- dependence of the interactions. 

We enumerate now the contributions to the strength functions Fr^d{kf) generated by 27r- 
exchange with virtual A-excitation. We start with (regularization dependent) contributions 
encoded in subtraction constants: 

FAk,)'-> = B,^, F.(M'"' = ||. (23) 

where the new parameter Bd = —M'^Vq{0)/4 stems from two-body Hartree diagrams and the 
momentum transfer dependence of the isoscalar central NN- amplitude Vc{q). Two-body Fock 
diagrams contribute only to F^-ikf) via a (subtracted) dispersion integral: 

Fr{kff'^^ = ^^ j^ rf/xIm(Vc + mc + 2/iVT + 6/i'W-T)|^ - ^ - QfJ^% 

+ ^ + 5/.^ arctan ^ + ^{2Ak) - ISk^ - ,^) In (l + ^) } • (24) 

The evaluation of the (left) three-body Hartree diagram in Fig. 2 leads to the results: 
,(3m_ 35g>| fl3 5 «^_35 / 5 3 _3\ ,, 



Somewhat more involved is the evaluation of the (right) Fock diagram in Fig. 2 for which we 
find: 

PrikfY''^^ = .,/-!,; , / dx 2Gs{x,u)Gs{x,u) +Gt{x,u)Gt{x,u) , (27) 

Gs{x, u) = {Qu* - 22u^ - 45 + SOx^ - 10m V) 

-|-4a;(10 -|- 9-u^ — 15a;^) arctan(-u + x) + arctan(-u — x) 

+ (35x2 ^ l^y^2^2 _ ^Q^4 _ 5 _ 9^2 _ 4^4^) j^ ^ ' ^^ ' -"^ ^ (28) 



1 + (m + a;)2 
1 + (u — xV ' 



u, ^,0/ 9x [1 + (m + a;)2l[l + (m — a;)^! 



4a; 16x2 

1 + C 

1 + (m — x) 



X (5 + 81^2 + ?,u* - l^u'x^ + 15x^) In | ^ {"" ^ ^^1' , (29) 
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.A,,.„.,3.,,„,:.4,.,.!<i-^l±^ ,30, 

with Gs,t{x, u) defined in eqs.(7,8). The strength functions F^,i{kf) are completed by adding to 
the terms in eqs. (23-30) the contributions from Ivr-exchange and iterated Ivr-exchange written 
down in eqs. (9,11,12,14,15,18,19,22,24,27,28) of ref.[Sn]- In order to be consistent with the 
calculation of the energy per particle E{kf) and the single-particle potential U{j), kf) we complete 
the Ivr-exchange contribution by its relativistic 1/M^-correction: 

FAkfY'""^ = ^"''f'' I ^u' - ^ + 2n(525 - 700u' - 96^^) arctan 2u 

-64m® + 744m^ - 1777^2 + M OSOm^ ' ^^ + -^j 1^(1 + 4^^) f • (31) 

The expression in eq.(19) multiplying the kinetic energy density T(r ) has the interpretation of 
a reciprocal density dependent effective nucleon mass: 



M*{p) =M 



2 



5k 

1 ^— + 2MFr(kf) 

28M2 ^^ ^' 



(32) 



We note that this effective nucleon mass M* (p) (entering the nuclear energy density functional) 
is conceptually different from the so-called "Landau"-mass M*{kf) defined in eq.(13). Only if 
the real single-particle potential has a simple quadratic dependence on the nucleon momentum, 
U{p,kf) = Uolkf) + p'^Ui{kf), do these two variants of effective nucleon mass agree with each 
other (modulo very small differences related to the relativistic (/c//2M)2-correction). 

In Fig. 7 we show the ratio effective-over-free nucleon mass M*{p)/M as a function of the 
nucleon density p = 2kj/3TT^. The dotted line corresponds to the result of ref.[36j based on 
Ivr- and iterated Ivr-exchange only. The full line includes in addition the effects from 27r- 
exchange with virtual A-excitation. It is clearly visible that the inclusion of the TrA^A-dynamics 
leads to substantial improvement of the effective nucleon mass M*{p) since now it decreases 
monotonically with the density. This behavior is of course a direct reflection of the improved 
momentum dependence of the real single-particle potential U{p, kf) (see Fig. 4). Our prediction 
for the effective nucleon mass at saturation density, M*(po) = 0.64M, is comparable to the 
typical value M*(po) — 0.7M of phenomenological Skyrme forces |SS1 IM]- The full curve in 
Fig. 7 displays another interesting and important feature, namely strong curvature effects at 
low densities p < 0.05 fm~^. They originate from the explicit presence of the small mass scale 
rriTr = 135 MeV in our calculation. 

Fig. 8 shows the strength function F^{kf) belonging to the (Vp)^-term in the nuclear energy 
density functional versus the nucleon density p = 2^^/37?^. The dotted line gives the result of 
ref.jHn] based on Ivr- and iterated Ivr-exchange only and the full line includes in addition the 
effects from 27r-exchange with virtual A-excitation. The subtraction constant Bd (representing 
density independent short-range contributions) has been set to zero. In the region around 
saturation density po — 0.16 fm~^ one observes a clear improvement since there the full line 
meets the band spanned by the three phenomenological Skyrme forces SIII j33]. Sly |S11 and 
MSk j32]. The strong rise of Fy{kf) towards low densities remains however. As explained in 
ref.jHE] this has to do with chiral singularities (of the form ■m~'^ and itl~^) in the contributions 
from Ivr-exchange and iterated Ivr-exchange. 
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Figure 7: The effective nucleoli mass M*{p) divided by the free nucleon mass M as a function of 
the nucleon density p. The dotted line shows the result of ref . [3^1 based on single and iterated 
pion-exchange only. The full line includes in addition the effects from 27r-exchange with virtual 
A-excitation. 
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Figure 8: The strength function F\/{kf) multiplying the (Vp)^-term in the nuclear energy density 
functional versus the nucleon density p = 2A;?/37r^. The dotted line shows the result of ref. [^2] 
based on single and iterated pion-exchange only. The full line includes in addition the effects 
from 27r-exchange with virtual A-excitation. 
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The knowledge of the strength function F^{kf) and the equation of state E{kf) allows one 
to calculate the surface energy of semi-infinite nuclear matter jSH] via: 



as = 2{?>Q'KpfY/^ rdp^pF^{kf)[E{kf) - Eo] . (33) 

Here, we have inserted into the formula for the surface energy a^ (see eq.(5.32) in ref-jSHj) that 
density profile p{z) which minimizes it. Numerical evaluation of eq.(33) gives a^ = 24.2 MeV. 
This number overestimates semi-empirical determinations of the surface energy, such as a^ = 
20.7MeV of ref.^U] or a^ = 18.2 MeV of ref.^^, by 17% or more.^ The reason for our high 
value Qs = 24.2 MeV is of course the strong rise of the strength function Fy{kf) at low densities. 
Its derivation is based on the density-matrix expansion of Negele and Vautherin [38^ which has 
been found to become inaccurate at low and non- uniform densities ^j. Therefore one should 
not trust the curves in Fig. 8 below p = 0.05 fm~^. Getting the right order of magnitude for 
Fy{kf) in the density region 0.1 fm~'^ < P < 0.2 fm~^ is already a satisfactory result. 

6 Nuclear matter at finite temperatures 

In this section, we discuss nuclear matter at finite temperatures T < 30 MeV. We are particularly 
interested in the first-order liquid-gas phase transition of isospin-symmetric nuclear matter 
and its associated critical point {pc, T^) . As outlined in ref . [H] a thermodynamically consistent 
extension of the present (perturbative) calculational scheme to finite temperatures is to relate 
it directly to the free energy per particle -F(p, T), whose natural thermodynamical variables are 
the nucleon density p and the temperature T. In that case the free energy density pF{p,T) 
of isospin-symmetric nuclear matter consists of a sum of convolution integrals over interaction 
kernels A^, multiplied by powers of the density of nucleon states in momentum space:^ 



dip,) - P^ 



27r2 



1 + exp 



p] - 2Mp 
2MT 



(34) 



The effective one-body "chemical potential" fi,{p,T) entering the Fermi-Dirac distribution in 
eq.(34) is determined by the relation to the particle density p = 4 J^ dpi pid{pi). We summarize 
now the additional interaction kernels arising from 27r-exchange with virtual A-excitation. The 
two-body kernels read: 

^(-) = 2A.W-^ + 20n^B,P-^ipl + pj) , (35) 



IT 



oo 



X:^^^ = - / dplm{Vc + 3Wc + 2p'Vt + Gp'Wt) 



2m^ 



P^+(Pl+P2)^ 4piP2 , 4]9iP2. 2 



X pin - r ^. -^^ + ^^(P^+P^) • (36) 



p2 + {pi - P2y p p3 

Note that the i?3-term in eq.(35) generates a temperature independent contribution to the free 
energy per particle, F{p, T)^^^^ = 37i'^B3p/2M'^, which drops linearly with density. Temperature 

■*One could reproduce the surface energy a^ = 20.7 MeV of ref.^^ by adjusting the short-range parameter 
Bd in eq.(32) to the value Bd = —75. The full curve for F^{kf) in Fig. 8 would then be shifted downward by 
29MeVfm^. Compared to B^ = —8 the fitted number Bd = —75 seems to be rather large. 

^ Since the temperature T is comparable to an average kinetic energy we count T of quadratic order in small 
momenta. 
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Figure 9: Pressure isotherms P{p,T) = p'^dF{p,T)/dp of isospin-symmetric nuclear matter at 
finite temperature T. The coexistence region of the hquid and the gas phase ends at the critical 
point: Pc ~ 0.053 fm-3, T^ ~ 15MeV. 



and density dependent Pauli blocking effects are incorporated in the three-body kernel /C3. The 
contributions of the Hartree and Fock diagrams in Fig. 2 to the three-body kernel read: 



4"' = ^^2p,p,(i + c) + 



2mlpiP2 



[ml + (pi + P2Y] [ml + (pi - p2f] 



m„ In 



ml + (pi +p2f 
ml + (pi -P2Y 



K 



(F) 



9\ 



4A/>3 



2X{p,)X{p2) + Y{p,)Y{p2 



X{pi) = 2pip3 - — - In 

2 ml + {pi - p3y 



(37) 
(38) 

(39) 



Yipi 



Pi 

4P3 



{5pl - 3ml - 3pl] 



I 3(pf - p^ + m^)^ + Amlpl ^^ ml + {pi + p^f 

ml + {pi- p^y 



IQpl 



(40) 



The remaining kernels building up the free nucleon gas part and the interaction contributions 
from Ivr-exchange and iterated Ivr-exchange have been written down in eqs. (4,5,6,7,11,12) of 
ref.jnj. It is needless to say that the extension of our nuclear matter calculation to finite tem- 
peratures T does not introduce any new adjustable parameter. 

Fig. 9 shows the calculated pressure isotherms P{p, T) = p'^dF{p, T)/dp of isospin-symmetric 
nuclear matter at six selected temperatures T = 0, 5, 10, 15, 20, 25MeV. As it should be these 
curves display a first-order liquid-gas phase transition similar to that of the van-der-Waals gas. 
The coexistence region between the liquid and the gas phase (which has to be determined 
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by the Maxwell construction) terminates at the critical temperature Tc. From there on the 
pressure isotherms P{p,T) grow monotonically with the nucleon density p. We find here a 
critical temperature of T^ ~ 15MeV and a critical density of pc — 0.053 fm~^ ~ Po/3. This 
critical temperature is close to the value Tc = (16.6 ± 0.9) MeV extracted in ref. |l2j from 
experimentally observed limiting temperatures in heavy ion collisions. In comparison, a critical 
temperature of Tc = (20 ± 3) MeV has been extracted in ref.jlHj from multi-fragmentation data 
in proton-on-gold collisions. Most other nuclear matter calculations find a critical temperature 
somewhat higher than our value, typically Tc ~ 18 MeV |^ 1^ 1^ . The reduction of Tc in 
comparison to Tc ~ 25.5 MeV obtained previously in ref. [H] results from the much improved 
momentum dependence of the real single-particle potential U{p, kfo) near the Fermi surface 
p = kfo (see Fig. 4). As a general rule the critical temperature Tc grows with the effective 
nucleon mass M*{kfQ) at the Fermi surface. 

The single-particle properties around the Fermi surface are decisive for the thermal exci- 
tations and therefore they crucially influence the low temperature behavior of nuclear matter. 
The inclusion of the chiral vrA^A-dynamics leads to a realistic value of the density of (thermally 
excitable) nucleon states at the Fermi surface. This is an important observation. 

7 Equation of state of pure neutron matter 

This section is devoted to the equation of state of pure neutron matter. In comparison to the 
calculation of isospin-symmetric nuclear matter in section 2 only the isospin factors of the 27r- 
exchange diagrams with virtual A-excitation change. We summarize now the contributions to 
the energy per particle En{kn) of pure neutron matter. The two-body terms read: 

^n(fc„)^^*^=i?n,3^+i?n,5^, (41) 

- 1 r°° ( Ak^ 

En{knT^^ = --I dpIm{Vc + Wc + 2p'VT + 2p'WT)\3pkr.--^ 

8vr^ Jam, I 3/i 

with i?„ 3 and Bn^5 two new subtraction constants. The relative weights of isoscalar (Vc^t) 
and isovector {Wc,t) NN-amplitudes have changed by a factor 3 in comparison to eq.(2). The 
three-body terms generated by the Hartree and Fock diagrams in Fig. 2 read: 



Enikn 



4 6 
.(3H) _ 9a^7t 



2 1 

-u^ + M^ - 3m^ + 5M^arctan2M - -(1 + 9m^) ln(l + 4^^) 



6A(27r/,)4 [3 4 



(43) 



Eniknf^^ = -{i^^, [d^ [G|(a:, u) + 2G|(x, u)] , (44) 

with Gs,t{x, u) defined in eqs.(7,8). We emphasize that in this section the meaning of u changes 
to M = A;„/m,r, where kn denotes the neutron Fermi momentum. All three three-body diagrams 
in Fig. 2 have now the same isospin factor 2/3 since between neutrons only the 27r''-exchange 
is possible. Note that the Pauli-exclusion principle forbids a three-neutron contact-interaction 
and therefore eq.(43) is free of any C-term. The additional contributions to En{kn) from the 
(relativistically improved) kinetic energy, from Ivr-exchange and from iterated Ivr-exchange have 
been written down in eqs. (32-37) of ref.j2|. 
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Figure 10: The energy per particle En{kn) of pure neutron matter as a function of the neutron 
density p„ = fc^/Svr^. The dashed curve gives the result of ref. P|. The full curve includes the 
TrA^A-dynamics and two adjusted short-range parameters Bn,3 = —0.94 and -B„,5 = —3.58. The 
dashed-dotted curve stems from the sophisticated many-body calculation of the Urbana group 
|16j. The dotted curve gives one half of the kinetic energy _Ekin(fc„)/2 = 3fc^/20M. 

Fig. 10 shows the energy per particle En{kn) of pure neutron matter as a function of the 
neutron density p„ = /cj^/Svr^. The dashed (concave) curve gives the result of our previous 
calculation in ref.[2l. The full curve includes the chiral TriVA-dynamics. The short-range pa- 
rameters Bn,3 and Bn.5 (controlling the contribution of a nn-contact interaction to En{kn)) 
have been adjusted to the values i?,„ 3 = —0.95 and i?„ 5 = —3.58.^ The dashed-dotted curve 
in Fig. 10 stems from the sophisticated many-body calculation of the Urbana group JH] ■, to be 
considered as representative of realistic neutron matter calculations. Moreover, the dotted curve 
gives one half of the kinetic energy of a free neutron gas, E]^i^{kn)/'2 = 3/c^/20M. Results of 
recent quantum Monte-Carlo calculation in ref.jlZj have demonstrated that the neutron matter 
equation of state at low neutron densities p„ < 0.05 fm~^ is well approximated by this simple 
form. One observes that up to p„ = 0.16 fm~^ our result for Enikn) is very close to that of 
the sophisticated many-body calculation jTHlEj. At higher densities we find a stiffer neutron 
matter equation of state. Again, one should not expect that our approach works at Fermi mo- 
menta larger than fc„ = 350 MeV (corresponding to p„ = 0.19 fm~^). One of the most important 
results of the present calculation is that the unrealistic downward bending of En{kn) (as shown 
by the dashed curve in Fig. 10) disappears after the inclusion of the chiral vrA^A-dynamics. This 
is a manifestation of improved isospin properties. 

^The short-range parameters Bn.3 and Bn,5 have been adjusted such that the asymmetry energy at saturation 
density takes on the value A{kfo) = 34 MeV. 
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8 Asymmetry energy 

As a further test of isospin properties we consider in this section the density dependent asym- 
metry energy A{kf). The asymmetry energy is generally defined by the expansion of the energy 
per particle of isospin-asymmetric nuclear matter (described by different proton and neutron 
Fermi momenta kp^n = ^/(l -F 5Y^^) around the symmetry line: 



EUK, A;„) = Eikj) + S^ A{kj) + 0(6' 



(45) 



Following the scheme in the previous sections we summarize the contributions to the asymmetry 
energy A{kf). The two-body terms read: 
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(47) 



In eq.(46) we have taken care of the fact that there are only two independent (S-wave) NN- 
contact couplings which can produce terms linear in density. It is surprising that also the 
other coefficient 10(3-B„^5 — -B5)/9 in front of the /c^/M^-term is completely fixed. This fact 
can be shown on the basis of the most general order-p^ NN-contact interaction written down 
in eq.(6) of ref. 48^. Out of the seven low-energy constants Ci, ... ,6*7 only two independent 
linear combinations, C2 and Ci + 3C3 + Cg, come into play for homogeneous and spin-saturated 
nuclear matter. The contribution of the three-body Hartree diagram in Fig. 2 to the asymmetry 
energy A{kf) has the following analytical form: 
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4 fi 2 
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- + AuA ln(l + 4u^) 
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(48) 



with the abbreviation u = kf/m.,^. The parameter (^ = —3/4 is again related to the three-nucleon 
contact interaction ~ {(g\/Af^) {NN)^ which has the interesting property that it contributes 
equally but with opposite sign to the energy per particle E{kf) and the asymmetry energy 
A{kf). Furthermore, both three-body Fock diagrams in Fig. 2 add up to give rise a contribution 
to the asymmetry A{kf) which can be represented as: 
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(49) 



The auxiliary functions Gs,t{x,u) have been defined in eqs.(7,8) and we have introduced a 
double-index notation for their partial derivatives: 



Gijk{x,u) 



x^u^ 



d^+^Gijx.u) 
dxWu'^ 



S,T, 



l<j + k<2. 



(50) 
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Figure 11: The asymmetry energy A{kf) as a function of the nucleon density p = 2^^/37?^. The 
dashed curve shows the result of ref . |2] . The full curve includes the chiral vrA^ A-dynamics. 



For notational simplicity we have omitting the arguments x and u in the integrand of eq.(49). 
The asymmetry energy A(kf) is completed by adding to the terms in eqs. (46-49) the contribu- 
tions from the (relativistically improved) kinetic energy, Ivr-exchange and iterated Ivr-exchange 
written down in eqs. (20-26) of ref. [2] . 

In the calculation of the asymmetry energy we use consistently the previously fixed short- 
distance parameters B^ = —7.99 and -B„.3 = —0.95, as well as -B5 = and Bn.5 = —3.58. 
Fig. 11 shows the asymmetry energy A{kf) as a function of the nucleon density p = 2k^r/3'n-^. 
The dashed (concave) curve corresponds to the result of ref.j2]. The full curve includes the 
chiral vrA^A-dynamics. The corresponding value of the asymmetry energy at saturation density 
po = 0.157 fm-3 is A{kfo) = 34.0 MeV. It decomposes as A{kfo) = (12.1 + 119.3 - 109.9 + 
12.5) MeV into contributions of second, third, fourth and fifth power of small momenta, again 
with a balance between large third and fourth order terms [2]. The value A{kfQ) = 34.0 MeV 
is consistent with most of the existing empirical determinations of the asymmetry energy. For 
example, a recent microscopic estimate in a relativistic mean-field model (constrained by some 
specific properties of certain nuclei) gave the value A{kfQ) = (34 ±2) MeV |inj. For comparison, 
other empirical values obtained from extensive fits of nuclide masses are A{kfo) = 36.8 MeV 
110] or A^kfo) = 33.2 MeV [30] • The slope of the asymmetry energy at saturation density, L = 
kfQA'{kfQ), is likewise an interesting quantity. As demonstrated in Fig. 11 of ref.jHIj the neutron 
skin thickness of ^^^Pb is linearly correlated with the slope parameter L. We extract from the 
full curve in Fig. 11 the value L = 90.8 MeV. This prediction is not far from L ~ 100 MeV 
quoted in ref. [UIj and L = 119.2 MeV obtained from the "standard" relativistic force NL3 
[S2]. Furthermore, we extract from the curvature of our asymmetry energy A{kf) at saturation 
density po the positive asymmetry compressibility K^^ = kJQA"{kfo) — 2L = 160.5 MeV. 

Again, the most important feature visible in Fig. 11 is that the inclusion of the chiral ttNA- 
dynamics eliminates the (unrealistic) downward bending of the asymmetry A{kf) at higher 
densities p > 0.2 fm~^ (as displayed by the dashed curve in Fig. 11). This is once more a 
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manifestation of improved isospin properties. 

9 Isovector single-particle potential 

In this section we generalize the calculation of the single-particle potential to isospin-asymmetric 
(homogeneous) nuclear matter. Any relative excess of neutrons over protons in the nuclear 
medium leads to a different "mean- field" for a proton and a neutron. This fact is expressed by 
the following decomposition of the (real) single-particle potential in isospin-asymmetric nuclear 
matter: 

U{p, kf) - Uj{p, kf) rs 6 + 0{6') , 6 = ^^^ , (51) 

Pn + Pp 

with U{p, kf) the isoscalar (real) single-particle potential discussed in section 3. The term linear 
in the isospin-asymmetry parameter 6 = {pn — Pp)/{pn + Pp) defines the (real) isovector single- 
particle potential f//(p, kf), and t^ -^ ±1 for a proton or a neutron. Without going into further 
technical details we summarize now the contributions to f//(p, kf). The two-body terms read: 

Uj{p, kf)^^'^ = (25^,3 - 53)^ + (2i?.,5 - B,)^,{k) + p^) , (52) 

k"^ r°° 

Ui{p, fc;)(2^) = -/- / dp lm{Vc -Wc + 2fi^VT - 2fi^WT) 



12vr3^2n^. 



The contribution of the three-body Hartree diagram in Fig. 2 and the three-body contact term 
can be written in analytical form: 

tV/(P,fc/j -9A(27rA)4^ l + (^-^)' [l + {u + x)^][l + {u-x)^] ^^^^ + ^^^/' 

(54) 
and the total contribution of both three-body Fock diagrams in Fig. 2 can be represented as: 
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(55) 



where x = p/mT^ and u = kf/rriT^. The auxiliary functions Gs,t{x,u) have been defined in 
eqs.(7,8). The (real) isovector single-particle potential Ui{p,kf) (restricted to the region below 
the Fermi surface p < kf) is completed by adding to the terms in eqs. (52-55) the contributions 
from Ivr-exchange and iterated Ivr-exchange written down in eqs. (28-36) of ref. [SB]- The imagi- 
nary isovector single-particle Wj{p, kf) (below the Fermi surface p < kf) has been discussed in 
section 4.2 of ref . [5^ . 

The generalization of the Hugenholtz- Van-Hove theorem J2L to isospin-asymmetric nuclear 
matter gives a model-independent relation for the isovector single-particle potential Ui{p, kf) at 
the Fermi surface {p = kf): 



Ui{kf,kf)=2A{kf) 
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Figure 12: The real isovector single-particle potential Ui(p,kf) as a function of the nucleon 
density p = 2/c^/37r^. The dashed and full curves correspond to the sections p = (bottom of 
the Fermi sea) and p = kf (at the Fermi surface), respectively. 



The second and third term on the right hand side just subtract non-interacting (kinetic) con- 
tributions from the asymmetry energy A{kf). We find at saturation density fc/o = 261.6 MeV 
an isovector single-particle potential of Uj{kfo,kfo) = 40.4 MeV. This is consistent with the 
value f/i ~ 40 MeV jHj deduced from nucleon-nucleus scattering in the framework of the optical 
model. The generalized Hugenholtz- Van-Hove theorem eq.(56) serves also as an excellent check 
on our analytical and numerical calculations. 

Fig. 12 shows the real isovector single-particle potential f//(p, kf) as a function of the nucleon 
density p = 2k^J3TT'^. The dashed and full curves correspond to the sections p = (bottom 
of the Fermi sea) and p = kf (at the Fermi surface), respectively. One observes a splitting 
of both curves which sets in at p ~ 0.10 fm~'^ and increases with the density. At saturation 
density po = 0.157 fm~^ the difference between the (real) isovector single-particle potential at 
p = kfo and p = is Ui{kfo, kfo) — Ui{0, kfo) = 7.1 MeV. This is much less than the analogous 



difference for the (real) isoscalar single-particle potential U{kfQ^ k 



fo) 



f/(0, k 



fo) 



26.4 MeV 



(see also Fig. 4). 

The full line in Fig. 13 shows the (real) isovector single-particle potential Ui{p, kfo) at satu- 
ration density kfo = 261.6 MeV as a function of the nucleon momentum p. In the region below 
the Fermi surface p < kfo the momentum dependence of this curve is weaker than that of the 
(real) isoscalar single-particle potential U{p, kfo) shown by the left half of the full line in Fig. 4. 
Finally, we note that the generalized Hugenholtz- Van- Hove theorem eq.(56) (which allows for 
an alternative determination of Ui{kf, kf) shown in Fig. 12) holds with high numerical precision 
in our calculation. 
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Figure 13: The full line shows the real isovector single-particle potential Ui{p, kfo) at saturation 
density kf^ = 261.6 MeV as a function of the nucleon momentum p. 

10 Summary and concluding remarks 

In this work we have extended a recent three-loop calculation of nuclear matter in chiral per- 
turbation theory by including the effects from two-pion exchange with single and double virtual 
A(1232)-isobar excitation. In the spirit of an effective field theory, we have encoded all short- 
distance contributions (from the high-momentum parts of the pion-loops integrals etc.) in a 
few adjustable contact-coupling constants. We have investigated a wide variety of nuclear prop- 
erties in this framework. The empirical saturation point of isospin-symmetric nuclear matter, 
Eq = —16 MeV, po = 0.157fm~^, can be well reproduced by adjusting the strength of a two- 
body term linear in density and weakening an emerging three-body term quadratic in density. 
The various constraints set by empirical nuclear matter properties (saturation point and poten- 
tial depth lead to the (minimal) parameter choice -B5 = and ( = —3/4, with little freedom 
for further variation. The momentum dependence of the real single-particle potential U{p, kf) 
is improved significantly by including the chiral vrA^A-dynamics. As a consequence the critical 
temperature of the liquid-gas phase transition gets lowered to the realistic value T^ ~ 15 MeV. 
The isospin properties of nuclear matter get also substantially improved by including the chiral 
vrA^A-dynamics. The energy per particle of pure neutron matter En{kn) and the asymmetry en- 
ergy A{kf) now show a monotonic growth with density. In the density regime p = 2p„ < 0.2 fm~^ 
relevant for conventional nuclear physics, we find good agreement with sophisticated many-body 
calculations and (semi)-empirical values. 

In passing we note that the inclusion of the chiral vrA^A-dynamics guarantees the spin- 
stability of nuclear matter ^^. These improvements can be traced back to repulsive two-body 
Fock terms as well as three-body terms with a very specific density and momentum dependence. 
Open questions concerning the role of yet higher orders in the small momentum expansion and 
its convergence remain and should be further explored. 

Our calculation takes seriously the fact that there exist two hadronic scales, the pion mass 
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rrin = 135 MeV and the delta-nucleon mass splitting A = 293 MeV, which are smaller than or 
comparable to the Fermi momentum kjQ ~ 262 MeV of equilibrated nuclear matter. Propagation 
effects of quasi-particles associated with these "light" scales are resolvable. Therefore pions and 
A-isobars must be included as explicit degrees of freedom in the nuclear many-body problem. 
Controlled by a systematic expansion in small scales [kj, 771-,^, A), the dynamics of the interacting 
TrA^A-system is worked out up to three-loop order. In this effective field theory approach the 
basic mechanism for nuclear binding and saturation are attractive 27r-exchange interactions 
of the van-der-Waals type on which Pauli-blocking acts in the nuclear medium. Most other 
phenomenological approaches ignore these " light" physical degrees of freedom and parameterize 
the relevant low-energy dynamics in terms of strongly coupled heavy scalar and vector bosons 
{a, uj, p, 6, etc.). Their propagation takes place on length scales of 0.5 fm or less and can 
therefore not be resolved in the domain relevant to nuclear physics. We are guided instead by 
a change of paradigm, namely that the nuclear many-body problem involves the separation of 
scales that is characteristic for low-energy QCD and its (chiral) symmetry breaking pattern. 

Appendix: Real single-particle potential above the Fermi 
surface 



In this appendix we summarize the continuation of the real single-particle potential U{p, kf) 
calculated in section 3 of ref. [Hj into the region above the Fermi surface p > kf. The expressions 
eqs.(8,9) in ref.jH] for the two-body potentials from the Ivr-exchange Fock diagram and the 
iterated Ivr-exchange Hartree diagram remain valid without any changes. Eq.(lO) in ref. [Sj for 
the two-body potential from the iterated Ivr-exchange Fock diagram gets replaced by: 
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(57) 



with the abbreviations u = kf/niTr and x = p/rriT^. Eq.(ll) in ref.|2| for the three-body potential 
from the iterated Ivr-exchange Hartree diagram gets modified to: 
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Finally, eq.(13) in ref.j3j for the three-body potential from the iterated Ivr-exchange Fock dia- 
gram is replaced by: 
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where we have again introduced the auxihary function: 

G{x) = u{l + u^ + x')-^[l + {u + xf] [l + {u- xf] In Yf'l% • (60) 

For the definition of the quantities ymin, s, a, t and R and the antisynimetrization prescription 
Ay we refer to section 4. 
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